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Abstract 

We introduce a method of rigorous analysis of the location and type of complex sin- 
gularities for nonlinear higher order PDEs as a function of the initial data. The method 
is applied to determine rigorously the asymptotic structure of singularities of the modified 
Harry-Dym equation 

Ht + Hy = -^H''+H''Hyyy I Ji" ( J/ , ) = ^ " ^ / ^ 

for small time at the boundaries of the sector of analyticity. Previous work 0, shows 
existence, uniqueness and Borel summability of solutions of general PDEs. It is shown that 
the solution to the above initial value problem is represented convergently by a series in a 
fractional power of t down to a small annular neighborhood of a singularity of the leading 
order equation. We deduce that the exact solution has a singularity nearby having, to 
leading order, the same type. 

1 Introduction 

The singularity structure of solutions of nonlinear partial differential equations in the complex 
plane is not a well understood subject. Insofar as the authors are aware, except for integrable 
cases, there are no general methods in the literature to analyze the type and location of singu- 
larities of solutions of nonlinear systems of PDEs in terms of the initial data. 

The goal of the present paper is to develop a relatively general and constructive technique to 
address this issue, which applies to equations that admit formal asymptotic solutions. 

In view of the complexity of the analysis, and for more clarity, we describe and apply the 
method on a number of concrete problems. It will be however transparent that the method is 
much more general. 

Formal asymptotic similarity solutions could exist for small or large time, or when one ap- 
proaches the finite blow-up time of a similarity solution of a PDE. We prove that complex sin- 
gularities of these formal asymptotic solution actually correspond to singularities of the solution 
of the full PDE. 
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The motivation for understanding complex singularity formation of PDEs, aside from intrinsic 
interest, is that in some cases of physical interest 0], there is evidence that singularities 
that appear in the real physical domain after a finite time can be traced to the complex plane. 

The procedure consists in the following steps: (i) an early time asymptotic expansion in 
powers of t, the validity of which is justified for the modified Harry-Dym equation in [Q, (ii) 
introduction of appropriately scaled "inner" dependent and independent variables beyond the 
region of validity of the expansion (i), (iii) determination of singularities of the leading order 
equation and (iv) proof that a secondary expansion in scaled time, involving inner- variables, is 
convergent in a domain encircling a singularity of the leading order solution. Insofar as analysis 
of the leading order equation (in step (iii) above) is concerned, which (typically a nonlinear 
ODE), formal calculations have been used before (see and references in JJ). These can now 
be rigorously derived from the general theory introduced in |14| . 

The present paper justifies the above four step procedure for the modified Harry-Dym equa- 
tion This equation arises in the small surface tension limit of Hele-Shaw interfacial evolution 
in the neighborhood of an initial zero of the derivative of an associated conformal map. The jus- 
tification of singularity formation is a crucial first-step to understanding "daughter" -singularity 
phenomena where a smoothly evolving interface corresponding to a zero-surface tension solution 
is singularly perturbed in 0(1) time by arbitrarily small surface tension effect. 

Consider the following initial value problem for the modified Harry-Dym equation: 

Theorem 36 and Corollary 37 in [S] imply that for any t e [0,r], for large enough \y — t\/t'^/^ 
with arg(?/ — t) € (—I"": f""); there exists a unique solution H{y,t) to (Q) with H{y,t) ^ y^^^^ 
as \y — t\/t^^^ oo, with the following asymptotic expansion for i <C 1: 

where = 1 and P„ is a homogeneous polynomial determined recursively in terms of Pn-i, 
Pn-2,---Pi- The first two polynomials are 

„ . ,^ 15 1, „ , ,x 25875 , 195 , 3,n 

Pi(a,6) = --a~-6, P^{a,b) ^ -^a" + —ah + -h" (3) 

Further, if we introduce the scaled variables 

'? = 7^ ; ^ = ^'"' ; i^(y(r/,t),t)=t^/'G(ry,r), (4) 

then, according to Corollary 37 in 5/, for \rj\ sufficiently large, with argr; G (^l"', I"")! the 
function G has a convergent series expansion in t: 

oo 

G(77,r)=^r'=Gfc(77) (5) 

fc=0 

^The variable variable C, = rf'l'^ is used there. 
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In this paper, among other results, it will be shown that the convergence of the series ^ 
actually holds in an extended domain in 77 that includes at least a region close to a singularity 
fjs of Go(?7) in a neighborhood of the boundary arg?7 g ("I""? I""") for large \fjg\. 

Substituting Q into we obtain the following equation for G{ri,T): 



G 2 

"9 -g'^^" 
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rGr 



— G'^ — G^G,pp^ — 



1 2 

gGo + gVG'o + GqG'q 



From ISJ, it follows that the leading order solution Go satisfies 



In order for G in jS)) to match the asymptotic expansion |(2Jl we need to require that 
Goiv) = ?7~^^^(1 + 0(1)); large, arg?7 G 



47r 47r 



(6) 



(7) 



(8) 



The solution Go(?7) to the leading order ODE ||7J) with asymptotic condition © have been studied 
before. Numerical solutions were found |2| and computational evidence suggested that there is 

a cluster of singularities rjs, where Go{r]) e"/^ i^)^^^ (^7 ~ Vs)"^^^- Using the fact that Go{r]) 
is indeed a similarity solution to the Harry-Dym equation, which is integrable, it was shown 
|S] that Q) can be transformed to Painleve Pn. Isomonodromic methods were used to prove 
existence and uniqueness of sectorially analytic solution for argry e ("I"''!''') that satisfies 
far-field condition (|HJ|. Outside this sector, the behavior of the transformed equation solution 
is given by elliptic functions, whose singularities can be related to the two-thirds singularity of 
Go (77), as above. 

However, unlike the isomonodromic method, the method based on generalized Borel sum- 
mation summation jl5| applies to initially small solutions of non- integrable equations as 
well. To apply this analysis in our example, which does not satisfy all the conditions in jl4|. 
small adaptations of the proofs are needed (see Appendix). One can determine that for large 
77, uniformly in the sector argyy £ [~ T" ~ T" ~ some 6 € (O, |7r), except for an 

exponentially small region around singularity of Go, 
form 

Go(ry) ^ r"^U{() 

where 



the asymptotic series of Ga{ri) is of the 
0{V~') (9) 



C = - log G + - log 77 



i4V2 



-n 



'9/4 + (2n - l)z7r 



(10) 



with the principal branch of the log, where G is a Stokes constant of Go in the large rj expansion 



for argTy g [—f'"' 



GoM - V-'^' 



1 + ^ Gkl] 



-9fc/2 



The function U (C) is determined implicitly from the equation: 



C = log4-2 



2Vc7-ln 



(11) 



(12) 
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and U{() has a singularity (where U = 0) at ( = (g = log 4 — 2 — itt, corresponding to a string 
of singularities at 77 = 77s, where 



iAV2 g/4 , 9 
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rii'" + - log ?7, = -2 + log 4 - 2hiTT + log C (13) 



where n g N has to be large for rjg to be large. For large jrysl (large h), it is to be noted that 
arg77s is close to — the anti-Stokes line. There is similarly another quasi-periodic array of 
singularities close to arg 77 = ^ , but our focus will be only on the ones in the lower- half plane. It 
can be shown that for large |7;s| the singularities of Gq lie within an exponentially small distance 
of Tjs and, to leading order, are of the same type. This can be further verified directly from the 
equation for Go- 

Remark 1 It is easy to check that Gq cannot be zero, except at a singularity rj — f]s. Furthermore 
, in any domain T) that excludes a neighborhood of the singularities of Gq, and extends to 00 so 
that arg 77 S [— f'"' + ^tt + 5'\, it follows from differentiability of the asymptotics of solutions of 
ODEs lEjl that 

sup|7;i/2Go(77)|, snp\r^y'G'i'm<C (14) 
Remark 2 By 7iear the singularity rj = rjs we have 

^[/3/2 + o(c/5/2) = C-C. = ^^i^l + '^vryv-Vs)[i + o{v^\v-Vsm (15) 

and hence for rj — rjs — o{rjs ^^^), for large enough \ris\, 

/ 5/6\ / „„. X 2/3 

Note that if 7'i|7;s|^^/^ < \ri — ris\ < ro\ris\~^^^ , with tq > small, then there exists upper 
and lower bounds for \U\, independent of rjs for large \ris\. Since the singularity f}s of Gq is 
exponentially close to rjg, it follows that the lower bound of Gq in this annular region is also 
independent of \r]s\- 



Given these leading order singularities for Go{ri), we investigate the series expansion (O 
known to converge for large enough \t]\ in any compact subset of arg 77 G (— |7r, |7r), in ; 
neighborhood of a singularity of Gq. 



Remark 3 The domain D in the next theorem, sketched in Fig. 1, is specified in Definition^ 
It contains a small annular region of a singularity rjs of U ( cf. f ) for large ti € N and 

a sectorial region arg 77 g [—^ + S,^ — S), for S small and large rj. 

Comparing powers of t'' (for fc > 1) obtained by substituting power series (O into results in: 

GlCkGk — Rk (17) 
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where the linear operator Ck is defined by 

and the terms Rk on the right side of (|17|) are given by 

Rk{ri) ^ - ^ Gk^Gk^Gk:,- ^ GkiGk^Gk^G'k^ (19) 

fci=fc-l kj<k,Y,kj=k 

In order to match to the asymptotic expansion expansion jSJ, we require 

Ak ( 47r 47r\ 

^fc(^) ~ ;^I+T72; 1^1 large, arg 77 6 (^- — ,— j (20) 

for some specific constants Ak (Ai = —1/2, A^ = |, A3 = — r^,...). As explained later, it is 
not necessary to impose (^0)1 : any solution Gk which approaches as I77I — > 00 with arg?7 e 
(— Itt — i5, |7r + (5) at a rate faster than 77^^/^ must necessarily have the asymptotic behavior 
PHl (See Remark El). 

Theorem 1 The expansion ^ is convergent in T> for all sufficiently small t . In particular, 
for any singularity fjs of Go{ri) near the anti- Stokes line argrj = — |7r with \ fis\ sufficiently large, 
there is a singularity of G{rj, t) for small r , to leading order of the same type, approaching it as 
r -> 0+ . 

Remark 4 The convergence of the Taylor expansion in t and the hounds on Gk and GJ, suffice 
to show that G{ri,T) has the singularities close to those of Go{ri) since for a circle S^-^ of radius 
ei around rjs we have 




For small t, G{r], r) thus has, to leading order, a branch-point of algebraic order 2/3. 
Remark 5 The convergence of the series (0) in T> is a corollary of the following lemma. 

> I, with A > I, < B < 1, so that 

(22) 
(23) 
(24) 



Lemma 2 There exist constants A and B independent of j 

\\V'^'G,\U^ < ^ 

] 

ll^j Woe,!) S — — 
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Remark 6 The proof of this key Lemma that leads to the proof of Theorem Q| is given at the 
end of § 6. First, we prove a Lemma bounding the Rk (rj) . This provides bounds of Gk using a 
suitable inversion of Ck in The estimates suffice for our purpose but are not sharp, as \2U\) 

implies a faster decay rate in rj. The uniqueness of the solution G{r],T) in the regime |r/| ^ 1 
for arg?] € (^—^n, |7r) with G{r], r) ^ rj"^/"^ is shown in 

The proof of Lemma[21is by induction; we first prove an general result for sums of type l|19(l . 

Lemma 3 With Gq the same as before, there is a constant so that for any A G (l,oo), 
i? G (0, 1), N 9 /c > 2 and {Gj}j=i^..._fc_i defined in T> which satisfy 1^2^) - ^^ we have in i^lUp . 

\\ri''/^Rk\\o..v < ^ {B^A^+BA"-') (25) 
Proof. It is convenient to break up Rk as: 

Rk = Ro,k + Rl,k 

where for k — I, 

Rom - -J- 

and for fc > 1, 

r> \ ^ i /^2y-^ o/^ r^'" \ ^ r^'" \ ^ 

-TtO.fc — -CzQ <-rfci Lf/C2 + 2"^0*-^fc-l ^ •^'^0*^0 "^fci(-rfc2 ~ ^0 / , tjfci (jfc2 "^fca 

+ 1 ^ .k 

" •^'^0 2^ '-^fci "-^fca ~ "J^O '^ki '^k2 <-^fc3 

where denotes summation over fc^ > 1 with — 

Ri,k{v) ^ 2 51 ^kiGk2Gk3 ~ Gk^Gk^Gk^G'lf^ 

The proof follows by using the upper bounds on Gj, G^ and Gj' in ^^-^^ for /c-1 > j > 1, 
using (|14|l and noting that 

— ^ k"^ ^ — T k'^ ^ — A k"^ ^ — A k"^ ^ — A k"^ I 

□ 
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2 Proofs 

The proofs rely on bounding Gk in (0. For given fco and 1 < fc < fco, it can be seen that the 
sohition Gfe to CkGk = that goes to as oo in the sector argrj e (— |7r + S, |7r — (5), 



with < (5 < is given by 



(26) 



Here 6*1 = — |7r + (5, 6*2 = fvr — (5 and 6^ — 0. Also, in (|26|l . ui, U2 and 1*3 are three independent 
solutions of CkU — 0, with the following asymptotic behavior for large (see 



exp 



exp 



4^/2 

4\/2 
^„ ? 
27 



,9/4 



9/4 



where is defined in l|18|l . (wi, f2, 1's)"^ is the third column of and 



$(?7) 



It is easily seen that for large |7y| we have 



Ui U2 U3 

u'l u'2 U3 
u'i u'i u'^ 



T] exp 



AV2 9/4 
27 ' 



V2iv) 



-5/8 



exp 



4\/2 

( 7 

27 



,9/4 



vsiv) ~ '''' ' 

and that wi, W2, W3 are three independent solutions to the adjoint third order linear equation 
C^v = 0, where the coefficients are regular when Gq 7^ 0. The Gk determined from H2()|) are 
bounded in any domain that excludes the singularities of Gq (the only places where Go = 0), and 
it is easily seen that the bounds H22|) - (|24|l are valid for 1 < fc < fco for large A and B (depending 
on fco). 

To prove the bounds (|22ll - 124l) in Lemma[21for all fc, it is sufficient to prove them for sufficiently 
large fc (large f3k)- 

Note. We have to treat separately two regimes of r/ and Pk ■ (a) |7;| = 0{el3^^^) or larger and 
(b) |7/| = o(/3^^^). These require different integral representation of Gk and choice of domain. 
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2.1 Control in the regime (a), |?7| > constant P^^^ 
It is convenient to define 

X = Pk^^^V and Zfc(x) = Gk{Pt^^x) 
Then, using variation of parameters (see §5), we have 

Zkix) = V [r] ix) 

where 

V[R]{X) = E 4 r - ^^(*)n,,3(x)i?(x)dx (27) 

where R depends on Rk and Zk] rij^^, and Wj are given functions of X; whose exact expression 
is irrelevant, with behavior n^^a = 0(x~^^^), W^i = " T X + o(l), W2 = - T ^'^ X + o(l) 
and W3 = 0(1) for large x; and P{{x): P2ix) and P3(x) are the three roots of the cubic 

«^ + ^x'/^a-X^/^ = (28) 
with the following asymptotic behavior for large x'- 

Pi = ^zx^/^-^lnx + o(l) 

P2 = -^^x'/'-llnx + o{l) 
^3 = ^In X + o(l) 

as X ^ 00 for argx S (-fi' + (5, fTi" - (5) , with < S < In ^7^, 61 = -|7r + (5, 6*2 = |7r - (5 
and 6*3 = 0. 

It is necessary that the operators V be defined in a suitable domain £ in the x-plane containing 
the integration path where the bounds for the previous Zj, j — 1, k — 1 are available to 
estimate Rk- Also V need to be bounded for large Pk- To satisfy the latter requirement for each 
j {j — 1, 2, 3), any point x £ ^ must have the property that it can be connected to ooe^^^ along 
a path Cj entirely in £ so that on the path x{s), parameterized by the arclength s increasing 
towards 00, 

^5RPj(x(s))> C\x{s)\'/' > 0, for j = l,2 

^ mm^)) > C |x(s)|-i > , 
as 

where C is a constant independent of x- It is shown in §3 that these properties are ensured if 
we choose 

X : X to the right of OSl, argx e (^-^'^ + - 
where OSl is the polygonal line connecting xi, X3 and X2, and where 

X3 = e , X2 - X3 + pe'2-/3 ^ xi^X3+ pe"^^/^ 
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Here p is chosen so that arg xi = ^ T" ^^"^ ^2 = ^ ^'^^ ^ is suitably small, independent 
of fc, so that 5 appearing in the proof of Theorem 1341 in §6 is smaller than i. The domain £ is 
sketched in Fig. 3. Corresponding to £, we define the domain £k (Figure 2) 



£k 



{■ 



/5-4/9 



xe£] 



2.2 Control in regime (b), rj = o (^P^ 
In this shown in §6, we can write 



,4/9 
k 



Gkiv) = v\Rk \ +^aj.gj(?7) 



where 



V 



5j = Go e'^^^k P ^ Rf. involves i?^, Gk and GJ, and 

rn I 

Piv) = / TT-r-^dil' , for some rji € V 



7?J ^^^^,Go(r7) / ■ Go(7y')i?fe('?')e"^''^"f^('')-^(''')ldV 



(29) 



(30) 



while uji — e*^^/"^, lu2 — ^ '^'^/'^ and 073 = 1 (the three cubic roots of unity). In H29|) the limits 
of integration satisfy rjj^k = Pk^'^Xji where xi, X2 and X3 are as defined in the last subsection. 
The choice of the domain Vk for -q is subject to the conditions below. 

1. Dk must contain a region So that winds around 77s, near argry = — |7r, excluding a 0(?7s ^^^) 
neighborhood of rjs itself. Since the singularity f}s of Go is within e~'-^l''=l distance of 77^ and 
a singularity of Gq is the only point where Go = 0, this condition ensures a lower bound 
on Go and provides the contour integration §g in Remark 2. 

2. Any point 77 G Dk can be connected to rjj^k along a contour Cj that lies entirely within Dk 
such that 5R[ajjP] is increasing monotonically when the points 77^,^ are approached. The 
integration contour Cj in 1)29(1 is chosen to be such a path. Monotonicity ensures there is 
no exponential growth in k {(3k)- A fortiori, the points rjj^k, as defined above, are points of 
maximum of [ujj P] in Vk ■ 

3. We must have for k > /cq, the property Vk+i U £k+i C Vk U £k- For k < ko, all Gk can be 
determined through the representation (|26|l on a common domain Vko U £ko ■ The necessity 
of this condition comes from the fact that the Rk , needed to determine Gk in the domain 
Vk i)£k, involve Gi, G2, Gk-i- 

4. For any k, Vk U £k must contain the region So around the singularity ris- Also, for large 
enough I77I in this domain, we must have arg 77 G [—f''' + (5, |7r — (5] . We shall furthermore 
ensure a nonempty common domain V = Hkyko V^k U £k]- 
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5. To simplify the analysis, we seek domains so that Dfc U £k does not contain turning points 
(occurring when arg 77 = ±|7r) of the WKB solutions for large k in igU), e^'=-f'^(x(''))+^^ (x('')) 
(see §4). 

3 Properties of P{r]) and choice of the domains V and Vk 

This section is devoted to the construction of the domains T) and Pfc corresponding to a particular 
77s, determined from H13|l for large enough n £ N. The domains clearly depend on the choice of 
77s. The construction is relatively involved since monotonicity of 3? [uJjP ] must be ensured, while 
P{rf) is only implicitly known through (p^- H12|l and (jSOJ- Also, such a domain has to contain an 
annular region around singularity ?7s. 

Remark 7 In this section, constants such as C , K , S, ri, r^, etc arc positive and independent 
of 7] and r/s • 

First we define T>A^ part of the region where Go{t]) ^ ri~^/'^. 

Definition 4 For small 6 E (O, ^) we have 

f , , 1, , 2 2 

'Dai = yr- M > i^lVsl, argr; e I --7r + <5, -TT- 



^ -{v ■ ^ = arg 77 e ( -^TT + (5, -^TT + (5 



9- . , 9- : ^^,\vim\>\v\>l\vs\, 

where 771 (6'(r)) — 2|7ys|e~*3'^+"' + T-e^'t for r > 0. We also define M = (— |7r + | and 

Va = Va, U Va, 

(See Figure 1). 

Lemma 5 For any point rj G T^a, there exist three piecewise smooth paths from 77 to 00, fj :— fjj , 
for j = 1,2,3 , contained in Va so that on any smooth segment we have 



d 



^iu;,P[ijis)])>Cm\'/'>0 



where s is the arclength. Furthermore, 

\7%s)\>cM>o 

Proof. On the line segment 7)(s) = 7;o + se*"^ in Va, H30|I and largeness of \fi{s)\ and \ris 
together with the asymptotic behavior Go (77) ^ fi~^/^, imply 



4- mcJ.P) = ^ 
ds 



vis) 



where 6j — argcjj £ {±|7r,0}, Ofj — arg77(s). For suitable and 770, it is easy to see that for 
any 77 G Va and j = 1, 2, 3, there exists a polygonal line so that cos (6j + ^9fj + 0) > C > 0. 
Further, the line can be chosen so that \fj{s)\ > \r]\. □ 
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Definition 6 



4 1 

L = { rj : argr] = --tt + S, -\t]s\ < hi < "^IVsl 



Definition 7 Let Sq be a region around r]s (the singularity ofU{({r])) in defined by 

So = [ri- r,|ry,|-5/4 < |r;-r?,| <ro|r?,|-5/4, arg(,7-r7,) G (_^ + ^,^+^)} (31) 

with < < ro, small enough to ensure argC/p^^ G [^|^' (^^^ relation flti\)). 
Definition 8 We define 

2?T,i = |^:k'/'-%^'l<3(^ii||J^ r for 770 G 5o and 

< Boc-^^l"''!''"'- < 1, for some r e [0, \ijs\^^'^] \ 



1 - Vu 



i + Vu 



For r G 



, we define 



(1) arg?7G j^argr/, - iiTslTy,! ^'^ + Ki\iq,\ ^'^r,~-T:^5, 



(3) |f7- 1| < 5Boe- 



-i<-4|r).L--ffi|r,,|^/-*r 



2(5 

2?T2 U ^T2,r, where J= Q—\r^,\^/^ 

Theorem 9 For large \ris\, for any point rj dV, there exist Bq, Ki and piecewise smooth paths 
from r/ to 00 fjj =: f] (j = 1, 2, 3^ contained in V, so that on any smooth subsegment we have 



ds 



^{uj,P[f,{s)])>C\f,{s)\"^>Q 



(32) 



Furthermore 



|77(s)| >Ci|r;| >0 (33) 
For the proof, given at the end of §2, we need a few more definitions, constructions and lemmas. 
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Remark 8 Note that by Lemma\^ it is enough to show that for any rj £ 'Dt, we can choose a 
path for each of j = 1,2,3 connecting to r]^ € L, entirely within 'Dt so that the monotonicity 
property is satisfied. Noting also that since the ratio of any two values ofrjG 'Dt is bounded 
by a constant independent of rjg, the second part of Theorem\^ follows. 



Definition 10 For k > ko, we define 



4/9 



Vi,k = VsM + Poe *3'','72,fc = ri3.k + Poe*^^, 773,^ = e/3\ 

where po is chosen so that argTyi ^ — — S, a.rgr]2,k — ^ + for < d < The parameter 
e is small, but independent of k, as needed in Lemma \25\ and fcp is chosen large enough so 
that for k > ko, we have e/?f^^ > M, for M as defined in Definition^ We define a boundary 
dEk = dE^ U dE^ where dE~^ is the straight line joining rj^^^k with 771,^ and dE^ is the straight 
line joining rj^^k to ??2,fc- We then define (See Fig. 2) 

Vk = V\£k 

Lemma 11 Given j = 1,2 or 3, for any rj G dE^, the path fj(s) from rj to rjj^k along dEk 
satisfies the monotonicity property ^32^) . 



Proof. On dE^ we note that 



—3? 
dp 



3? 



e'27r/3y^l/2 



for some positive constant C. On dEi^ we note that 



dp 



e-*2V3^i/2 



|77|i/2sin 



TT 1 



|l/2 



for some positive constant C. It is therefore clear that the path 773(3) from 77 to 773 ^ satisfies 
the monotonicity property (|32|l . In a similar manner, it is seen that 3? [wi^*] and 3? [ci^2^'] satisfy 
1)321) on a path from 77 to rjj^k for j — 1,2 along dEk- □ 

The lemma above, together with Theorem^prove the following Corollary: 

Corollary 12 (Property 1:) For all sufficiently large k, given any point rj e 7)^, there exists 
a piecewise smooth path Cj for each j = 1,2,3 from rj to r]j,k such that the path is entirely in Vk 



(34) 



^5ft(a;,P(77(s))>q77|i/2>0 
Furthermore, if fj £ Cj, we then have \ fj\ > C\ri\. 

Remark 9 To prove Theorem\^ we introduce three autonomous flows as follows. 
Definition 13 Let rij[t,rjo) be the solution to the differential equation 



r] = ie-'"t'^Lo-^Go{il), where lui = e*2'^/^ 



-i27r/3 



W3 = 1 



(35) 
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with initial condition 77^ (0,770) = ?70j where (f>j are given by: 



TT Gtt 2 

<Pl = 4>2 = —, 4>-i = gTT (36) 

Remark 10 We note from I^SO\} that for any choice (pj E (0,7r), 

^ 5ft [u;,Pirj,it,rjo))] = cos (0, - |) > 
Hence using arclength parameterization we have 

when fi[s) G T) . Thus, the differential equation generates ascent paths for ^[LOjP]. 
Lemma 14 There exists a Bq so that So C Vt- 

Proof. Since for 77 e Sq, the corresponding [/ determined from (|12|l has upper and lower bounds 
independent of 7]s, as discussed in Remark|21 Also, from for 77 £ So, aigU^/^ e [— f''', f t!"] • 
Thus, it follows that for 77 g 5o, we have |1 — \/[7|/|l + Vu\ < Bq for some Bo < 1. Thus, for 
some Bo, we have Sq C I?t,i C X't- □ 

Definition 15 It is convenient to define, see IS5\) and 1,^6]) . 

Vj = ie-"^^uj-^ 



Remark 11 It follows that 

TT 137r TT 

argz/i = --, argz^2 = — , argi/g = (37) 

The specific choice of (and thus of Vj) is unimportant, but it is essential that <j)j, &TgVj remain 
in compact subintervals o/(0,7r) and (— |7r, -l) respectively, independent of rjg and 6. 

In order to study the solution to (|35|l near 77^, it is convenient to think of U{t) ~ U{r]{t)) as 
an unknown together with 77(t). Using (|10|l . (|12|) and H35|) . it follows that 

2 1^3/2 ^ [1 ^ ^^(^)] ^^(^) ^ ' ^ (38) 



[/ = -aj\r^s?'^Vu{U - 1) [1 + £^1] [1 + E2] , where 

£^2(r;) 



i^/2 



aj|77s 



3/4 



,3/4 



8^3/2 



(39) 



13 



where Ei, E2 will be shown to be small for large \'qs\ in the range of integration and 



The initial condition Uq satisfies 
«4y2 9/4 



i^/2 



,3/4 



3/2 



27 



27 



- 2x/C/( 







(40) 



Remark 12 /i is 6e noted that with (f)j given by Hclb]} and using the fact that as n ^ 00 (i.e. 
o-s \ris \ — > 00 j, we get argry^ ^ It follows that in this limit, 



3' 



argai - 

It is important for us that argaj £ §) 



IOtt 
'21 



arga2 — -, argaa 



(41) 



Lemma 16 For suitable and 5 < if Vofl ^ I^T^y then for sufficiently large \rjs\ and some 



te (^0,^\ris\^^'^^, 77j(i;?7o,o), leaves Vt^ through L. 



Proof. The differential equation satisfied by 77 and the corresponding ?7 for < t < 2-^\r]s\'^/'^ 
is given by: 



,3/4 



(;7-i)(i + S4) 



where 



Es^ EiU+{U -1), Ei^ {VU -l) + VUEi, 
where Ei, E2 are defined in (|38|l and (|39|) . It follows that 

v'^'^Vo!o+l'^j [\^ + E3)dt 



{U-l) = {Uo^o - 1) exp 



iV2 



,3/4 



8r/3/2 



{1 + E4)dt 



(42) 



(43) 



where J7o,o is obtained from H4U|I by substituting 770 = ?7o,o- It is convenient to define the leading 
order equations 



f^^^Vo,o+2''^t ; f/-l = (C/o,o-l)exp 



^ -3/4 



8773/2 



dt 



From (gSl and (gijl it follows that 



,3/2 _ -3/2 ^ 



E^dt 



(44) 



(45) 
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C/-C/=(C/o,o-l)< exp 







1 exp 









1^2 



3/4 



8773/2 



exp 



(1 + E4,)dt 



3/4 



8773/2 



(46) 



From iB^ it follows that 

argryo,o + Stlry^l"^/^ > argry > argr/0,0 + 7;Kzt\-qs\^^''^ where 2X3 min sin ( ^tt + i/j) 

Using these relations in H44|) we have 

|(C/-l)| = |(C/o,o-l)|e-^^^l''=l'''* 
where 2Ki is a lower bound (independent of 5) of 



2|r?,|3/4 



iV2. 



,3/4 



8773/2 



for 77 restricted to the domain \fj\ > ^\ris\, arg77 G [— ^tt — (5, — |7r + 4(5] . Thus, for some t in 
<t < 2^1^8 1"^^^ J ''7 leaves the domain Vt^ through the segment of L, when 5||77s| < Irjl < ^\r]s\. 
Now, we show that 77 is close to 77 and hence has roughly the same behavior. We define 



(C,v^) = 



7^3/2 _ ^3/2 



Uo,0 - 1 



on the interval [0, ii], for < ii < ■^It?^^/^. We introduce the norm 

|l(C,y)||oo- sup m\+ sup ei^'^M^'^'lVit)] 

0<t<ti 0<t<ii 

and consider the right side of 145() and 14611 as the mapping 
of the ball 

B = {m,V{t)):\\{C,V)\\^<e,} 

for some small ci in the Banach space of pair of continuous functions {C{t),V{t)) of t in the 
interval [0, ii] for ti < ^Ivsl^^'^- 

Using the smallness of E3 and E4 for large 77 it can be checked directly that 
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and that 

II (^i(Cl,yi),^2(Cl,^l)) - (^l(C2,1^2),^2(C2,V2)) lU < e2||(Cl,^l) " (C2,X^2)|U 

for some €2 < 1 and the map is contractive. Thus, there is a unique solution to the integral 
system l)45|l - (|46|) for {({t),V{t)) in B. In particular, this implies that 

Pit) -1\< |C/o,o - l|e-^^l"=l'^**, Mt))'/' - 4I^,\ < 3t (47) 
Hence, with r as in the definition of I'T2,r we have 

|C/-1| <5Boe-^^l''=le-^^l"»l'''(*+'-) 
argr; > argr;o,o + Kslrjsl-'^^h > arg?/, - ifshsl"^/* + K^lvsl'^^/^it + r) 
|^3/2| g |^|^^|3/2 _ K,\rjs\'/' - 3(t + r), \rj,\'^' + K,\rjsp' + 3(< + r)) 

Therefore, from the definition of I?t,2, for small enough t + r, we have 77 € T^t,2, while from 
continuity, there exists some larger t + r < -^|?7sP^^ for which 77 S L as it exits 25^2 • n 

Lemma 17 Let 770,0 G 2?t,i- Define 

Then, there exist Bq and Ki so that 1) G Dt^ U I?t2 for large \rjs\ and < t < |?7s|^^^- 
Proof. Note that for any Q <t< \r]s\^^'^-, we write (EU and ((23) as 

rf"^ = + \f^ ^jU{l + E^)dt , = 60,0 exp ^-a,\r^s?'^ j\l + i?2)(l + 

^ ° ^^^^ 

where |6o^o| < -806"^"'''°' ^ , with chosen in accordance to LemmalT^and 2Ki := min^ cos (-1 - 

We introduce Tj{t) and U{t) (describing leading behavior) by 



i + Vu 

It is to be noted that 



i = bo^oe-^M^''^ , f'' = 4^ + %, Uit')dt' 



^3/2 _ -3/2 = ^ r [{U - 1) + f/£;i]dt (49) 
2 Jo 



\/C/ 1 - V C/ , -a t\v \^'* 

^ = 6o,oe "^''"'l 



i + %/c7 1 



exp|-a,|7y,|3/4^ [(1 + £;2)(1 + ^i) - l]rfi| - 1 



(50) 



We note that is a lower bound for 5R [ctj] for I775I large. It is convenient to define the pair 
of continuous functions. 



16 



and the norm 

||(C,y)||oo- sup \m+ snp ei^^l^^l'^'^lFWI 

0<t<ti 0<t<ti 

forti G (0, |77,|i/''). Consider the right hand side of and as a mapping (T^ (C. V).T',(C. V)) 
on the baU 

B^{{C,V) : ||(C,V)||oo <eiti} 

Using smaUness of Ei, E2 and their derivatives with respect to 77, it can be readily checked that 
(^1,^2) is a contractive mapping of the ball B into itself; hence the solution (C, ^) satisfying 
H49II and H5U|I is in B for large \rjg\. In particular, since SRa^ > |if4 we have 

< Soe--^l-l^^^(*+^), Mt)f'' - ^lT\ < 3 + r) 

There are two cases: if t + r < |77s|^/'', then clearly 77 G I'ti- If ^ t + r < 2\ris\^/'^, from 

the definition of Vt^, it follows e X't,2, with = 6 (j^) • □ 

Proof of Theorem 1^1 From Lemmas 151 II 71 [T^ (see Remark |S1 as well), it is clear that that the 
domain V = Vt U Va is invariant under the flows fjj{s). From Remark ITUl Theorem |51 follows. 



1- VW) 



4 Properties of Pj{x) and choice of domain E 

Remark 13 The WKB solution for large f3k of the homogenous equation CkU = (see 
item 5) is not uniformly valid in the domain T) for large r\. To invert the operator Ck in the 
regime rj = 0(/3^^^), we introduce the scaled variables: 

X = (51) 
The WKB solution to the homogeneous equation is then of the form 

QhP,{x)+w,(x) ^here a = P' are roots of the cubic + ^ax^^^ - X'''^ = (52) 

9 

We now choose a domain £ where the WKB solution is valid. First, we define a boundary DEl, 
which corresponds in the x plane to dE^ (see Definition \10\) . 

Definition 18 Let OEl = |x • V = Pk^^X ^ dE}^. We define d£^ and d£^ analogously in 
terms of dE'^ and dE^ , (see Definition ] lUi) . 



Definition 19 We let 

X : X to the right of dE^ , arg r/ e 
(See Fig. 3.) It is also convenient to define 



27r , 27r ^ 

h(5, S 

9 ' 9 
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Remark 14 Note that for large k we have the following properties : T) C X'fe U £k and 2?fc+i U 
Ek+i C Pfc U ffc. This follows from the construction of Vk and £k- Our strategy is to prove the 
bounds in Lemma\^in the domain Vk U£k based on bounds on all previous Gj, j = l,2...,{k — 1) 
established on the domains Vj U £j (which contain 2?^ U £k)- The large k requirement is not 
restrictive, since for any fixed kg it is possible to choose A large enough so that the bounds in 
Lemma\^ hold for 1 < j < feg. 

The main theorem in this section is the foUowing. 



Theorem 20 For any x ^ ^ j possible to choose a path Cj connecting x to ooe ' , where 

f + 5, 02 = 



?i — + d, 62 = — 6 and 63 = so that, except for a finite set of points. 



and 



where s is the arc-length increasing towards 00 and the (different) constants C above are inde- 
pendent of X- Furthermore, for \x\ sufficiently large in £, and with x G Cj as above, we have 
\x\ > C \x\ for C > independent of x, and 

Proof. This follows, after a few Lemmas, at the end of §3. □ 

Remark 15 Though the domain £ restricts the size of\x\ (it is bounded below), it is convenient 
to first consider the properties of Pj on an enlarged domain £q with no restriction on \x\ and 
larger width: 



Definition 21 



£0= <x- argxe 



9 ' 9 



It is convenient to associate each Pj with a first order differential equation as follows. Note 
from ^ that with Fj x^^^V" 

we have 

X-^/^ = V' + ^V (53) 
Now, we consider the trajectory in the complex x plane generated by the differential equation 

The solution with initial value xo will be denoted by Xi(^; Xo)- Using H53|l . it follows that 

^-(2 + 9^-2)2 



dt 4(2 + 27i/;2) 



(55) 
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For large x G £ it is clear from 1)531) that the three possible behaviors of V' are ip ~ *\/ 1' 

ip = —i\J\ and '0 ~ |x~^^^- We associate these behaviors with P[, and Pg respectively, so 
that 

P[ - * {Ix"" ,P'.- - ^ {{x"" , - \x-' (56) 

Remark 16 Note that ds — \^\dt, and so on a trajectory generated by the differential equation 
i54}) , we have 

^ 5RP,(X) = |P,'(x)l 

and hence one of the two conditions in Theorem\2U^is satisfied by the pathCj — {x ■ X — Xj 
provided it remains within £ . 

Lemma 22 5RPi increases monotonically on the boundary of £q counterclockwise from ooe^s^ 
to ooe~*9^ with 

-f 5RPi(x(s)) > C\x{s)\^'\ 
ds 

while 5RP2 increases monotonically on the boundary of Eq clockwise from ooe^^^'^ to ooe^^'^ with 

^ mixis)) > C|x(s)|^/^ 
ds 

s being arc-length on Sq. 

Proof. Consider the solution to H55|l . with initial condition on the imaginary 0-axis slightly 

above = "^yf^- This corresponds to starting at x = cx^e*^^/^ with P{(x) and tracing the Stokes 
line where Im Pi = and 5RPi is increasing. From the equation it is clear that remains on 
the imaginary axis and approaches ioo, implying that argx = |7r is a Stokes line where SftPi is 
increasing monotonically all the way to the origin in the x-plane. This also means that locally 
near x = 0, P' ^ i^iX"^''^ and Pi |tt'ix'^^^, since this is the only root of the cubic (|5^ which 
is real on x = re'^'^^^. This corresponds to ■0 '^iX~^^^ as x ~^ 0. Now, taking the initial 
condition slightly above — i'lj^, it is clear from the differential equation (|55|l that ip remains 

on the positive imaginary 0-axis and approaches i/j = i from below. This corresponds to the 

fact that arg x — — |7r is a Stokes line beyond the turning point X — Xs — i j e^'^'^' ^ , with 
5iPi increasing monotonically towards ooe^'^^^ and for large r, ^ JRPi > Cr^^*. Now, consider 
the segment x = ^ g-427r/9^ where < r < ( 4"^^ J • If we introduce the transformations 



into lO, then 

2 

- + q^^ = 0, where q = r"^!^ 
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The roots of the cubic that corresponds to Pj are: 

* - = -^''''ITj'^W - . rn-. where J^l- .[7^^ (57) 

^ 18gi/3 1 3(2916)1/3 ji/3 J V 59049 ^ ^ 

(the principal branch is used). The asymptotic behavior of ^'j for small r is given by 

From ||57I), it follows that on the hne x = re"'^'"/^, for < r < ( ) , we have 

d 



dr 



iV2 J 

3?Pi(re-*2'^/9) = r^/^SR^-i > Cr^f^ > 



Thus, for all r, we have ^5RPi (re '27r/9j ^ Cr^/^. From the reflection-symmetry between Pi 
and P2 on the positive real ^-axis, the statement for P2 follows. □ 

Lemma 23 diP^ decreases monotonically on the boundary of 80 counter-clockwise from c«e**9^ 
to 0, and 



ds |x(s)P/2 + l 

s being the arc-length towards 00. In this, the positive real x-a2;is is a Stokes line with SftPs 
increasing towards ooe**' and satisfying the above monotonicity condition. 



Proof. Consider (|55|l starting with on the positive imaginary axis, slightly below -0 = iy 

corresponding to x = ( ) e~*^'^'^. The differential equation implies that ip remains on 

the positive imaginary axis as it moves towards the origin. This corresponds to % = ooe~'o'^, 

since ip ~ fx^^^'' for large Xi where P3 ~ Thus, the segment x = re"*^'^, r > i^^^j is 
a Stokes line with 

^nPs (re-^%A > - 
dr \ / r 

From the symmetry about the real x-axis, the same argument can be repeated for x = re*^'^ 

/ /oN4/9 

for r > ( 1 to show that this segment is also part of the Stokes line with 5RP3 increasing 
with r. 

For r < ( ] I an examination of in (|57|l shows that 5RP3(re** a'^) continues to 
decrease monotonically with decreasing r, though these segments are not part of any Stokes line. 
Near the origin, given the asymptotics of in 158II . it follows that P3(x) ~ fx"^^^- Hence a 
corresponding inequality follows, incorporating this behavior at the origin, while at the same 
time satisfying condition for large x 

which implies the inequality in the Lemma. Furthermore, when argx = 0, it is easily seen that 
Pg is real and positive and so P3 increases monotonically to 00 as we approach ooe*°. □ 
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Lemma 24 For any S € (0, ^) there exists Rq independent of S so that 
1. 

as 

for C > independent of any parameter, where s is the arc-length representation of part 
of the boundary of £ for which |x(s)| > Roi s is increasing in r when x — re~^9'^+^^ and 
decreasing when x = re'a'^"*''. 

2. 

^5iP2(x(s)) > C\x{s)\'/' 

for C > independent of any parameter, where s is the arc-length representation of part 
of the boundary of £ for which \x{s)\ > Roj s is increasing in r when x — re's'^"*'' and 
decreasing when x = re"*^'^^"'. 

3. 

^m{xis))>C\x{s)\-' 
as 

for C > independent of any parameter, where s is the arc-length representation of part of 
the boundary of £ for which \x{s) \ > Rq; s is increasing in r when x = f exp {±z [|7r — i5] } . 

Proof. This follows from the asymptotic behavior of P[, P2 and Pg for large x in H56|) after 
noting that 

□ 

Lemma 25 For < ei < r < Pq- There exists a small enough 5 > 0, independent of any 
parameter, so that 

-^5RPj (re"'i'^+'*) > C> for j = 1,3 



5RP2 (re-^i'^+'M > C > 



while 

_d_ 

dr 

with C independent of 6. Again, for ei < r < Rq, there is a S > 0, independent of any parameter 
so that 

d 



— 5RP, (re'i'^-*'^') > C > for j = 2, 3 
while 

-^5RPi (re*§^-**) > C> 

for some C independent of 5. 

Proof. From the lemmas about the behavior of P, on d£, the statements are clearly true for 
6 = 0. From continuity, it follows that the same is true (adjusting C) for all sufhciently small S 
and hence the lemma follows. □ 
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Definition 26 

where 



OSl = dSt U SZ 

dSl = {x = X3 + re'^'^/^ for < r < |x2 - Xsl } 

{x = X3 + re-*2"/3 for < r < |x2 - Xi } 



Lemma 27 SRPs increases in r on d£^ and dSj^ . JtPi decreases in r on d£^ , but increases in 
r on d£7 ■ 3?-P2 increases in r on d£t and decreases in r on d£7 and in all cases, we have on 



d£ 



L, 



5RP,(x(r)) 



> C > 



where C only depends on the choice of\xj\. ^Pj attains a maximum on d£L at the corresponding 
Xj- 

Proof. We note that since jxsl is small, we have 

-^mixir)) = -5? [J'^(x(r-))e'^'^/^] ~ |x(r)|^/^sin J + > C> 
where arg x = 6' G + ~ By symmetry we also get for x on d£2 



dr 



mixir)) = -SR N(x(r-))e-^^'^/^l ~ IxWI^/' sin ( J - ^ ) > C> 



6 2 



For Pi we find that for x € d£^, 



d 



-£SRPi(x(r)) ~ \xir)\'/' cos (- + -)> C> 



dr 



3 2 



On d£j^ , we obtain 



dr 



5RPi(x(r))~|x(r)|^/'cos 



>C>0 



Thus, on OSl, ^Pi increases monotonically from top to bottom with Jj3?Pi(x('S)) > C > 0. On 
this boundary P2 increases monotonically from bottom to top by a similar argument. On the 
other hand, P3 is maximum at Xs; it decreases as we move up or down. □ 

Lemma 28 On the boundary of£, KPi increases monotonically with s as we traverse the bound- 
ary counterclockwise and: 

^KPi(x(s))>C|x(s)|^/">0 
whereas JtPj increases monotonically with the arclength s as this boundary is traversed clockwise 



^3?P2(x(s)) > C\x{s)\'/' > 
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On the other hand at the upper part of d£ , i.e. on d {£ f^ {x '. Ini x > 0}), 



ds 



5RP3(x(s))>C|x(s)ri >0 



where the boundary is traversed counterclockwise. For the lower part of d£ , i.e. ond {£ f^ {x '. Im x < 0}) 
we have 



d_ 

ds 



mix{s))>C\xis)\-'>0 



where the boundary is now traversed clockwise. 
Proof. The proof follows from Lemmas I54I57I □ 

Proof of Theorem 1201 Any x ^ 9£ can be joined to ooe*^^ along d£ so that ^3?Pj(x(s)) 
satisfies the lower bounds given in Lemma If x <= ^, we choose steepest ascent paths for 3tPj 
until (i) it goes to oo, or (ii) it intersects d£, from which point we continue along the ascent 
paths of d£. The proof is complete. 



5 Estimates on the solution Gk in the domain Sk 

The main theorem proved in this section is the following. 
Theorem 29 For rj G £k we have 



QO,£k — 1 



felloe 



< K 



Gq 
Rk 



Gq 



oo,£k 



ca,Sk 



\\Gl\ 



< K 



Rk 



where the constant K is independent of k (and therefore of [3]^). 

Remark 17 The proof comes at the end of §5, after a few lemmas. It is convenient to derive 
an integral equation for Gk and its first two derivatives. We exploit the largeness of k to control 
the solution of the integral equation. The asymptotic behavior of the solution of the homogeneous 
equation CkU — is different in the regimes rj fc**/^ and \ri\ > fc**/^. Consequently, different 
integral equations will be used in £k and Vk (analyzed in the next section). 



In £k, it is convenient to introduce scaled variables: 



Then, H17() becomes 

^kZk 



'^-fe+/3fc 



*0 , *1 



X 



/3feX^ 



Zk 



^0 
Gf' 



kpA/i 

3 ' 




(59) 



(60) 
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where 



and and are defined by 



1 



3/2 



^0 SGqGq" 

~~3" 



77 



*1 



(61) 



(62) 



From the large ij G £k behavior of Go we see that and ^'i are bounded for large Pk as well 
as for large x- Let v be the solution for x € f of 



ChV = TZ 



(63) 



Using rigorous WKB results ^Hli it follows that for large (3k, there exist three independent 
solutions of the associated homogeneous equation, with leading behavior vi, V2, fa where 



where a = Pj are the three roots of the cubic equation 

Note that two roots of coincide iff + wX^^^ = i.e. iff 



(64) 



(65) 



X^Xs 



( 81V3\ 



4/9 



3±i27r/g 



only possible outside £. Hence the Vi,i = 1,2,3 are independent in £. The corresponding Wj 
are given by 

iP'P- 

(66) 



W'- 



J 



and the Pj are uniquely determined by the following asymptotic conditions for large x'- 

Pi = ^ix'" - ^ Inx + o{l),P, = -^^X''' - ^ Inx + o(l), P3 - ^ Inx + o(l) (67) 



Wi 
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lnx + o(l), W2 
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lnx + o(l), W^^o{l) 



(68) 



We now use the Vi to write an integral equation for w, equivalent to 1)63(1 . with appropriate decay 
conditions at 00. First, we have 



M 



Vi V2 V3 



where for j = 1, 2, 3 



1 1 1 

77121 '7122 "723 

77I31 77132 777,33 

W' 



^2 

7;3 



(69) 
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From the asymptotic properties of Pj and Wj , it follows that for large Pk we have rriij 
for all Furthermore, for large x, we also have 

^21 = 0(X'/'),"J22 = 77123 = 0{x-') 

msi = 0(x'/'), r7i32 = 0(x'/2^, 77733 = Oix'^) 
Qi = {M' - Q2M)M-\ where Q2 = Pk 



Let 






1 











1 




_2^5/2 






Then A4 satisfies the differential equation 

M' -iQ2+Qi)M = 

Denoting 



nil ni2 ni3 

7721 ?^22 ^23 
"31 «32 "■33 



1 1 1 

77721 W22 77723 
TO31 777-32 77733 



and 

we have 



A = r7722m33 - r7l23"^32 - m2im33 + r7l2im32 - 7773177722 + m3im23 



"1,3 = ('"23 - m22)/A; 772,3 = (^21 - 77723)/A; 773,3 = ("^21 - m22)/A 



The first two rows of — Q2M. are zero. Hence, the same is true for the first two rows c 
Therefore, 

'0 " 
Qi=Pk^ 

^31 ^32 &33 



Defining rj so that 
we have 



rj = SP^wf + 3P;W^' + SW^P^' + P^" + (wf + SW^W^' + W^") 
We note that Vj =0(1) for large /3fe. For large x we have 

ri-0(x-'/'), r2 = 0(x-'/'), rg = 0(x-') 

Also, with Ai = (77721 — 77722)(77).32 — 77733) — (77722 — 77723)(77731 — 77732)we have 
^3,2 = [{ri - r2)(77732 - 77X33) - (?'2 - r3)(77l3i - 77732)]/Ai 
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^3,3 = -[(''I - r2){m22 " m23) - (r2 - r'3)(m2i - m22)]/Ai (81) 



h,! = 7-3 - &3,2"l2,3 - &3,3"^3,3 (82) 

From the asymptotics of and m^.j for large /J^ we get b^ j — 0(1). For large x G ^ we have 

63,2 = 0(X"'), ^3,3 = 0(x"'), and 63,1 - 0(x"') (83) 
Then, for large x ^ follows that 

ni,3 = 0(x-'/'), n2,3 = 0(x"'/') and 713,3 = 0{x^^^^) (84) 

In the domain £ the b^ j are analytic in x, bounded for large /3fc and decay for large Xi (see (I83||V 
Fm'thermore, H73|l implies 
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Vi := Wj- —v^ 



f3k ' 



gPkX - 03,2 



Using variation of parameters, we see that one solution of CkV — TZ satisfies: 

vix)^V [r] ix) ; v'ix) = V' [r] ix) ; v"ix) = V" [r] ix) 

where 

i?(x) - nix) - - h,2v'{x) - Pkb^Mx), 

and the operators V, V' and V" are defined by: 

3 



[P.(x)-P.(x)]+w,(x) - ^^(^'>n,,3{x)R{x)dx 



Pfc Jo. 



3 



nj,i,{x)R{x)dx 



V'mx) = Vm3,,(x) / [^.(x)-P.(x)l+l^.(x) - 3(x)i?(x)dx 



(85) 

(86) 
(87) 

(88) 
(89) 
(90) 



where 9i 



jTT — (5 and 03 — 0, and the paths of integration Cj are chosen 



to be the ascent paths for ^Pj of Theorem EOl Also, note that for large Xj W^j(x) grows at 
most logarithmically with x implying that, uniformly in £, we have Wj = o{PkPj)- As we shall 
see, there is a unique solution v of (|SS|l that decays as x ^ 00 in with R having similar 
decay properties. The way we show this is by proving contractivity of the integral system in a 
suitable space of decaying functions. (In fact, there can be no other decaying solutions, since 
the associated homogeneous equation does not have nonzero decaying solutions in £.) 
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Lemma 30 // the Pj satisfy Property 1 in £ , then for sufficiently large k ( or, which amounts 
to the same, large f3k ) we have 



Pk 



c 



c 



oo,e < ll-R||oo 
Pk 

C 



l!V"[i?]|loo,£< \\R\\oo 
Pk 

where the constant C is independent of R and Pk ■ 

Proof. Theorem 1201 shows that on Cj (defined before Lemma On|) we have |x| > C\x\ and 

^5RPi,2(x(5)) > C\x\'/\ ^^^P^im) > C\x{s)\-' 
Since / P' is bounded, this imphes that for sufliciently large we have 



(91) 

(92) 
(93) 



d 
ds 

d_ 

ds 



Pi, 



f3k 
Pk 



Also, fi-om ||7nj) and (|7T)l . 

|m2,i| < C\x\'/\\m2,2\ < C|x|'/^ \m,.s\ < C \x\-\ 

M < C |m3,2| < C |x|^/^ |m3,3| < C 1x1"^ 

while from In^J < C\x\^-'/^. Then, 



PI 



00c 3 



exp [PkiPAx) - PAx)] n3^,P(x)e'^^(>^)-^^(*)dx 



< 



C\\R\\oo\x\ 
PI 



-3/2 /■! 



d[exp(/3fe[5RP,(x)-5RP,-(x)])] (94) 



The bounds for V follow; V' and V" are bounded similarly. □ 

Corollary 31 Define the operator Tk acting on triples (zk, z'/^, z'^) as follows: 



Tk {zk,z'k, zl) (x) = -Pl^h^,zzl - + 63,2 ) 4 + h 



X^ PkX' 



h,! Zk 



(95) 



Then, it follows 

llx^/'v [Tk {zk, 4, 4)] \U.£ < c IPk^klUs + Pk'\\x'^'4\U.£ + Pk'Wx'^'zkWoox 
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II V" [n {zk,zl4)] < C [p^'KlUe+P^'Wx'^'zlW + ||x'/'^fc||oo,£ 

Proof. This follows from Lcmmar^and bounds on 63^ in (|83|l and those on ^fo, 5*1 that follow 
from jnH). □ 



Lemma 32 



x'^'v' 

V" 



< 



GO,£ 



< 



GO,£ 



PI 

PI 



Rk 



Rk 



Gq 



00, f 



00 ,f 



00, s Pf^ 



Rk 



Gg 



Proof. This is a consequence of LemmaOOl with R replaced by 0l/^ Rk/GQ. □ 



Lemma 33 For H^^Gq'^II^^ < 00, and for Pk sufficiently large, the system 1^86]) has a unique 
solution {zk{x)j z' (x), z" (x)) in £ , which satisfies the hounds 



Pk 



,4/3 



Rk 



G^ 



00, £ 



11^5/2^; II c. < ^ 

llA -2fclloo,£ :i „2 
Pk 



4/3 



Rk 



GO,£ 



ll4'l|oo,£ < 



CPk 



A/3 



Pk 



Rk 



Gq 



OD.S 



Proof. Define the Banach space !F of triples of functions {u,v,w) analytic in the interior of £ 
and continuous in its closure in the norm 

\\{u,V,w)y = /3,'/'||x'/'"l|oo,£ +/3,'/'||x''/'«l|oo,£ +/9,:'/'|k||oo,£ 

We associate Zk, z'f, and z^' with u, v and w respectively, and consider R as depending on u, v 
and w for fixed Rk/Gg. We define the linear operator L : T ^ J- hy 

L[{u,v,w)] := (y R{u,v,w) , V' R{u,v,w) , V" R{u,v,w) 

where V, V', V" are now thought of as acting on u = z^, v = z'j^ w = z'^ for fixed R^G^^. From 
Ht)U|) . (|85|) and (|87|l . the definition of 7^ in l|95|l . and the estimates in Corollary 1311 it is easily 
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seen that 

\\L[{u,v,w)]-L[{u,v,w)]\\jr 

= II V [Tk{u ~ u,v — v,w — w)] ,V' [Tk{u — u, V ~ V , w — w)] , V" [Tk{u — u,v — v,w — w)] 

C_ 
PI 



Hence L is contractive and the system (|86f) has a unique solution (z^;, z^, z'^). The estimates on 
z/j, z'f,, z'l foUow easily from Lemma l32l □ 

Proof of Theorem 1291 This is a consequence of Lemma noting that 



6 Estimate of Gk for large k in 

In this section we prove the following. 
Theorem 34 In I?^ (see Definition MU^) we have 



\v'^'Gk{v)\\oo,v, < 



Rk 



Gq 



oc,'DMEk 



Rk 



Gl 



(96) 
(97) 



Remark 18 The proof of theorem {34\ ) is completed at the end of ^6, after a few lemmas estab- 
lishing the properties of C^^ . 



We first find a representation of the solution to 



7k- 1 



Rk 



C-kGk — Gfc' + -TrpyjVG'k „ „Q Gk + —pp-Gk — t^t 



9G: 



9Gg 



Go 



(98) 



for large fc for 77 G where |?7| is small compared to Again following P!^, there exist three 

independent solutions ui, M2, U3 to the homogeneous equation LkU = such that, for large 
we have 

Ujiv) ~ gjiv) = Go{-q)e'^=^'>'"'P^'^\ where wi = e'2'^/^ uj2 = e"''^'"''^, loj, = 1 and 



An 

We use gi, 53 to find a suitable integral equation for the solution u to H98|) . As in §5, it is 
convenient to define 



M := 



9i 92 93 

(ik"'9[ p-k"'9'. Pk'^'9's 



pf^'9'{ pr''9'i nr-'g'i 



ij-iii 0-2/3 , 1 



(100) 
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and 



h := {M' - Q2M)M-\ where Q2 := 13, 



1/3 






1 3G^' 



1 

1 

^ 



We get 

Using we see that 



9GJ /3feGo 9f}l'^Gl 

M' -[Q2 + Qi)M = {) 



-2/3 






O ^o'-'O I /I 'jQ 






2 

"G^ Go ^ 9GJ 



1/3 



The cohimns of Ai also satisfy (I102|l : it follows that for j=l,2,3 we have 



^k9j ■■= g'j 
We note that 
where 





2GI 




Go 



53,2 



9', + 1 


-4- 

1^ G% 


= Cku 






277 


\ Go 






Go 



Pfc ^0 '-^0 . ^^0^0 

) LfO Lrg 



TQ^-ro 



For large |?7| in 7)^ we find 



(101) 



(102) 



(103) 



(104) 
(105) 



&3.2 = 0{if/^), 634 = 0(77"3) (106) 

Also, 634 and 63,2 are analytic in Vk- It follows that the Gk in (|98|1 also satisfy the integral 
equation 



where 



Gfe(ry)=v[i?,] {fl) +Y.a,g,{7^) 



Rkivi) = 7^3 - ^3,2Gfc - 63.lGfe 
Go 



(107) 
(108) 



The constants aj are defined in piU|) in terms of Gk{i]i^k), Gk{ri2.k) and Gkiflz.k) and the 
operator V is defined by 



V[RkM = Y.^^^Go{v) / Go(r/)i?.(ry')e-^''^''[^('')-^(''')ldry' 



(109) 
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The contours of integration chosen in (|109|l are ascent paths of 9? [ujjP] , see Coroharv ll2ll . Given 
Gkim,k), Gk{ri2,k) and Gkim.k) we define ai, 02, 03 by 



Again, 



where 



and 



giivi.k) g2{m.k) gsim.k) 
9i{.V2,k) 92{m,k) gaim.k) 
gi{m,k) g2{v3,k) gaiva.k) 





ai 






0.2 






. "3 . 





Gk{m,k)-V[R][7^i,k] 

Gk{m,k)-^[R][m,k] 

Gfe(r/3,fc)- V[i?][ry3,fc] 



V [i?] [ri) - 5] ^c.,7i,(ry)Go(77) / Go{n')R{ri') 



hj{ri) 



Gq 



It is to be noted that 



\r^-^'^hM\<C 



for some constant G independent of 

A few properties of V and V' follow from Property 1 of P{r]) (established in §3). 



(110) 



(111) 



(112) 



(113) 



Lemma 35 Assume WRWoo.Vk < Then, 



h^/^V [R] (^)IU < ^||i?|U 
Pk 

for a constant Ki independent of (3k ■ 

Proof. Note that on any of the contours Cj, from Property 1, there exists a constant G > so 
that fj > C\-q\ for i) G Cj and 

j-^^{to,P{^{s))}>CM^)\"^>^ 
where s is the arc length. Therefore the proof follows from the estimate 



Pf^'G,{^) / Go(ry')i?('7')e''^''("^[^(''^-^('''^"rfry' 

1 



< 



G 



l!^l|oo,-Dfc 



□ 
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Lemma 36 Assume |l-R|loo,X'fc < oo. Then 



where K2 is a constant independent of /3fe . 



Proof. As before, there exists a constant C > so that on the contour Cj we have C\rj\ < rj' 
and 

d 



-5ic.,P(7y'(,s))>Ch'(s)|i/2>o 



ds 



where s is the arc length. Thus 



< / d{exp[m{p{v)^piv'm 

Jo 

The Lemma follows by noting that in Vk we have I??] < |?73.fc|. □ 

Corollary 37 We have 

\\rj'/'V[bM + bs,iG,] < ^ \\\v'/'G',\Uv, + W^'G^lUv, 

Pk L 

Proof. This follows from Lemma 1^ and the bounds on 63^2 and 63^1 in (|106|l . □ 
Corollary 38 We have 



ll^lloo.Pfc 



\v'^'G',\\^^v, + \\v'/'Gk\U 



Proof. This follows from Lemma 15^ and the bounds on 63,2 and 631 in (|106|l . □ 



Corollary 39 The following inequality holds 

'Rk 



< 



Rk 



Gl 



ca,Vk 



for a constant independent of k. 
Proof. This follows from Lemma 1551 □ 



Corollary 40 We have 



Rk 



< 



K, 



3/2 



Rk 
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Proof. This follows from Lemma 15^ after noting that for rj G Vk, \ri\ < "q^.k- □ 
Definition 41 Define the linear operators Ti and Ti by 

3 



3 



(see 1^1 10\) ) since r]j,k € d£k, Gkirjj^k) are known from the previous section. 



Lemma 42 FKe have 



||r?3/2Ti[G,,G'fe]|U,p. < 



2{K + Ki) 



Rk 



Gq 



oo,-DfcU£t 



felloo 



Proof. From (|110|l . since rjj^k are large and therefore gjivj' ,k)/gj{ilj.k) are exponentially small 
in f3k for j' j, it is clear that 



— I — fc3,2Gjj — bz,iGk 



From Lemma H35|l. Corollaries 1371 and 1381 it follows that 
Wjv'!k9jiv,,k)\ < 2\ij,,kf^'\Gk{Vj,k)\ + 



2Ki 


Rk 






Pk 


Gl 




f3k L 



Now, we conclude from Theorem 1291 that 



{Vom) - Gk{V],k) 



\v'^'Gk\Uv, + \\v'^'G',\Uv 

(114) 



h,k\'^'\Gk{v,.k)\ < f 



OQ,£k 



(115) 



Since rff'^gj{ri)/{rij^^gj{7]j^k)) are bounded independently of j3k and the proof follows. □ 
Lemma 43 



W^%[Gk,G'^]\\oo.v, < 



Rk 



Gq 



oo.VkUEk Pk 



Proof. Taking into account the behavior of hj{rj) for large rj we note that 

\Pk^^V^j!khj {rij,k)ajgj {Tjj,k)\ < G Pl'^\iij,k?/^\7i/^^ a^g^ (77^- fe)| 
Using (|114|l and (|115|l . the proof follows. □ 
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Proof of Theorem 1341 We consider the space Banach B of pairs of analytic functions {u, v) 
in the interior of Vk continuous in its closure with the norm 



\iu,v)\\ 



1773/2^11 



■,■0,- 



Associating Gk and GJ. in H1U7|) and with u and v, we define the linear operator L from B 

to B by 



L[{u,v)] - V Rk[u,v] +Ti[u,v]y 



Rk \u,v\ 



where Rk is now thought of as an operator on (u, v) for fixed ^ such that Rk[Gk, (77) equals 
the right hand side of H1U8|) . 

It is a simple application of Lemmas 1351361 142I43I and Corollaries 1371 and 1401 that 



\\L[iu,v) - iu,m < m^^y) - iu,v)\\ 



where 



max ■ 



3/2 r> 3/2 ■ 
K3 Ks K4V3,k KsV3.k 



1 

<2 



f3k' Pk' (3l/' ' pl^' 

for sufficiently large (3k and small e. Contractivity of L implies that it has a unique fixed point. 
The estimates in the Lemma follow from Hl()7|l and 

Proof of Lemma 121 First, for k — 1, ko, the statement in the Lemma holds if A is sufficiently 
large (depending on fco) in a common domain VkoUSka i chosen to contain Vko+i^^ko+i- Assume 
therefore that fc > fco where fco + 1 is large enough to ensure contractivity in Theorems 1291 and 1341 
Assume the statement holds j — 1, fco in a common domain Vk^ U £ko a-^d for j = fcg + 1, 
fc — 1 in a corresponding sequence of domains Vj U £j . It follows from the construction of these 
domains that it then holds in Vk U £k- We then get the estimates on Rk needed in Theorems 1291 
and 1341 which imply 



k\\oD,VkU£, 



< 



K 



10 



Ri 



h'/'G^lloo,I5.U£. <i^ll 

c 

\W'^Gk\Uv.U8, < 



Rk 



r<3 



and therefore, from the estimates on ||_Rfe77'^/^|| in (|25|l . we get 

K10K3 



i„5/2r" 



k\\oo,VkU£k 



< 



fc3 
K11K3 
fc2 



{B^A'' + BA''^^) 
{B^A'' + BA''-^) 



Using eq. H17|l and the bounds on Rk, it follows that 

\\Gl'\Uv,u£, < K^,\\r^^/^Gk\Uv,u£, + K,3\\ri'^'G[\Uv,u£, 



K 



13 



K3 



14 



fc2 



{B^A" + BA" 



It is clear that for B sufficiently small and A sufficiently large, the estimates (I22f) - (I24II on Gk, 
G'f. and G'^' follow. The result follows now by induction. 

Proof of Theorem ^ Now this follows easily from Lemma 12 since the estimates guarantee 
convergence of the Taylor series (O for sufficiently small t. 
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7 Appendix: Singularities of nonlinear ODEs 



We first mention briefly a number of results in |14| and then allow for slight modifications in the 
assumptions, to adjust for the equation of Gq. 

7.1 Setting of jl4| and generalizations 

We adopt, with few exceptions that we mention, the same conditions, notations and terminology 
as |15) and jl4| : the results on formal solutions and their generalized Borel summability are also 
taken from "15". 

The differential system considered has the form 

y' = f(x-\y) yeC", xeC (116) 

where 

(i) f is analytic in a neighborhood Vx x Vy of (0, 0), imder the genericity conditions that: 

(ii) the eigenvalues Xj of the matrix A = — |^^(0, 0)| are linearly independent over 

Z (in particular A^- 7^ 0) and such that argAj are all different. 

We now allow for the same assumptions, except we replace (ii) by 
(ii') There is at most one zero eigenvalue of A and all the other Xj are linearly independent over 
Z (in particular Xj 7^ 0) and such that argAj are all different. 

By elementary changes of variables, the system (|116() can be brought to the normalized form 

y'--Ay + -iy + g(a;-\y) (117) 

X 

where A — diag{Aj}, A — diagjaj} are constant matrices, g is analytic at (0, 0) and g(a;~^, y) = 
0{x^'^) + 0(|yp) as x 00 and y ^ 0. 

As in we normalize the system so that ^{ctj) > 0. 

Performing a further transformation of the type y i— > y — X^fcLi ^kx^^ (which takes out M 
terms of the formal asymptotic series solutions of the equation) , makes 

g(N-i,y)-0(a;-^^-i;|y|2;|x-2y|) (x ^ ex.; y -> 0) (118) 

where 

M > max3?(aj ) 
j 

and 0(a; 6; c) means (at most) of the order of the largest among a, 6, c. 

Our analysis applies to solutions y(a;) such that y(a;) — > as a; ^ cxd along some arbitrary 
direction d = {a; G C : arg(a;) = 4>}. A movable singularity of y(a::) is a point a; G C with 
x~^ G Vx where y(a;) is not analytic. The point at infinity is an irregular singular point of rank 
1; it is a fixed singular point of the system since, after the substitution x — z~~^ the r.h.s of the 
transformed system, ^ = — z~-^f(z,y) has, under the given assumptions, a pole at z = 0. 
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An n-parameter formal solution of pi7ll (under the assumptions mentioned) as a combination 
of powers and exponentials is found in the form 

y(a;) = ^ C^e-^-^"" x"-^si,{x) (119) 

ke(Nu{o})" 

where Sk are (usually factorially divergent) formal power series: §0 = yo and in general 

00 

^>^(-) = E ^ (120) 

that can be determined by formal substitution of (|119|l in (|117|l ; C € C" is a vector of parame- 
ters^(we use the notations C'' = 11^=1 ^ ^ — (^Ij ■■■T^n), ot = (ai, ...,Q!„), |k| — fci + ... + fc„). 

Note the structure of (|119|l : an infinite sum of (generically) divergent series multiplying 
exponentials. They are called formal exponential power series 

From the point of view of correspondence of these formal solutions to actual solutions it was 
recognized that not all expansions H119(l should be considered meaningful; also they are defined 
relative to a sector (or a direction). 

Given a direction d in the complex x-plane the transseries (on d), introduced by Ecalle |10| . 
are, in our context, those exponential series (|119|l which are formally asymptotic on d, i.e. the 
terms C^e~^'^''x°''^x~'' (with k e (NU {0})", r e NU {0}) form a weU ordered set with respect 
to ^ on d (see also QHI)-'^ (For example, this is the case when the terms of the formal expansion 
become (much) smaller when k becomes larger.) 

We recall that the antistokes lines of IjllTI) are the 2n directions of the x-plane iXj R+, —iXj K+, j = 
l,...,7i, i.e. the directions along which some exponential e~^^^ of the general formal solution 
H119|l is purely oscillatory. 

In the context of differential systems with an irregular singular point, asymptoticity should 
be (generically) discussed relative to a direction towards the singular point; in fact, under the 
present assumptions (of non-degeneracy) asymptoticity can be defined on sectors. 

Let d be a direction in the x-plane which is not an antistokes line. The solutions y{x) of 
H117() which satisfy 

y{x) ^0 {xed; \x\ 00) (121) 

are analytic for large x in a sector containing d, between two neighboring antistokes lines and 
have the same asymptotic series 



y(2:)-yo {xed; ^ 00) (122) 
In the context of (|117|l . a generalized Borel summation CB of transseries H119|l is defined in 

m- 

■^In the general case when some assumptions made here do not hold, the general formal solution may addi- 
tionally logs iterated exponentials, and powers 1101 . The present paper only discusses equations in the setting 
explained at the beginning of the present section. 

^We note here a slight difference between our transseries and those of Ecalle, in that we are allowing complex 
constants. 
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The formal solutions (|119|l are determined by the equation pi7|) that they satisfy, except for 
the parameters C. Then a correspondence between actual and formal solutions of the equation 
is an association between solutions and constants C. This is done using a generalized Borel 
summation CB. 

The operator CB constructed in jTQ can be applied to any transseries solution H119|l of (|117|l 
(valid on its open sector S trans, assumed non-empty) on any direction d C S trans and yields 
an actual solution y = CBy of H117() . analytic in a domain San- Conversely, any solution y(a;) 
satisfying 1122|l on a direction d is represented as CBy{x)^ on d, for some unique y{x): 

y{x) = J2 C^e-^ '^-x^ '^y^ix) - ^ C^e-^ '^^x^ '^CBykix) = CBy{x) (123) 

k>0 k>0 

for some constants C G C", where Mj = [^oij\ +1 ([-J is the integer part), and 

oo 

y'^(-) = E = M) (124) 

(for technical reasons the Borel summation procedure is applied to the series 

- x''°'sk(:r) (125) 

rather than to Sk(a;) cf. ifTH^jl . ifT^ ^ . 

The modification necessary to extend H123|l to the case Aq = is outlined in tl7.3l 



7.2 Normal form of Eq. (|7D 

We first give some detail on the normalization procedure, in the limit |a:;| oo. It can be checked 
that there is a one-parameter family of formal solutions to Q in the form Cx^^^^ — x~^ + .... 
The physical problem requires C = 1; this suggests the substitution Go = x^^^^ + h{x) where h 
is expected to behave like — -^a;"^. 

The normalizing substitution produces an equation with solutions in the form (|123|l . where 
the terms with k > contain exponentials with argument linear in the final variable; the type 
of the exponenential in the equation for h can be found by linear perturbation theory around a 
solution ho] with h — Hq = S, the leading order equation for 5 is 

S'" + ^x'^^H' + ^x^^H = (126) 

where the substitution of the form S — A{x)e'^^'' shows that p — 9/4 implying that the natural 
variable is x^/^. 

Taking Go = x~^/^ + x^^/^g{x^/*), C = x^/* in we obtain 



1 „ / 11 32 1 
which, written as a system, becomes 









/ 243 ' 







0(9=, r") (128) 




32 
729 













1 
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The eigenvalues of the first matrix on the rhs of (|128|l . {0, ±^^^}, are the values of A in H123|) . 
The fact that one eigenvalue is zero requires a slight modification in the proofs of |15j . 

7.3 Extension of the proofs in |15| to the assumption (ii') 

In an attempt to minimize the possibility of confusion with the setting in I15j we assume that 
the order of the system is n + 1, we count dimensions starting with zero, and take Aq = 0. There 
is no contribution from Aq to the general formally decreasing transseries 1)1191) ; this is due to the 
normalization ^{aj) > 0. 

The convolution equations satisfied by Y = By and Yk = Byi^ are given still given by 
equations (1.13 ) and (1.16) as in Jl] (with the notation A = —B used there): 

-pY = Fo-AY-i?PY+7V(Y) (129) 



(-P + A - k • a) Yk + (b + k • m) T'Yk + ^ dj * Y*j = Tk 

|j|=i 

(130) 

The only difference relevant to |2| with respect to the analysis in ^Sj is in the study of 
Y, and once the analog results are obtained, the analysis of Yk is virtually identical. By the 
normalization choice, we have Fq — p^^H(p) where H is analytic at zero. In the equation (2.35) 
of UHl 

Y=(A-p) ^ (Fo-BPY+AA(Y)) =M(Y) (131) 

we separate the zeroth component which is apparently singular (as was done in the study of Yi 
in|15j §2.2.2; here the analysis is simpler): 

- p(Y)o - ao7'(Y)o = Fq^o + := Rq (132) 

or 

-p(Y)[,-(Y)o-ao(Y)o=R^ (133) 
which we rewrite as an integral equation, which after integration by parts reads: 

(Y)o = -Fo;o + (1 + «o) Fo;o(ip)dt - (AA(Y))o + (1 + ao) {M{Y))^ {tp)dt = A^W(YX134) 

Jo Jo 

The system is of the form ()131|l 
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Y = X[il(Y) (135) 

with A^[^l = A4 for all components other than the zeroth one defined in 11341) . The equation 
p35|) is contractive in the ball B = {Y : {p : \p\ < e} : ||Y||oo < 2e} for small enough e, and also 
in the focusing algebra (3a) in §2.1.1 in for /3fe = 1 (allowed by the normalization of Fo) as 
follows from immediate estimates. 

No other nontrivial adaptations are needed in the proofs in |15| . 

7.4 Results of fT^ as extended in ^7.31 

The map y i-^ CB{y) depends on the direction d, and (typically) is discontinuous at the finitely 
many Stokes lines, see Theorem 4. 

For linear equations only the directions \j E_|., j = 1, n are Stokes lines, but for nonlinear 
equations there are also other Stokes lines, recognized first by Ecalle. CB is only discontinuous 
because of the jump discontinuity of the vector of "constants" C across Stokes directions (Stokes' 
phenomenon); between Stokes lines CB does not vary with d. 

The function series in (|123|l is uniformly convergent and the functions yt are analytic on 
domains San (for some (5 > 0, i? = R{y{x),5) > 0. 

Theorem 44 There exists 6i > so that for |<^| < di the power series 

oo 

F™(^) = '^ = 0,1,2,... (136) 

k=0 

converge. Furthermore 



{^{x)) {xeSs,, x^(x) (137) 

m=0 

uniformly in Ss^ , and the asymptotic representation J j5'7| ) is differentiable. 

The functions are uniquely defined by | |i,y7| ), the requirement of analyticity at ^ ^ 0, and 
ro(0) = ei. 

Remark 19 A direct calculation shows that the functions F,„ are solutions of the system of 
equations 

^Fo -r' (AFo-g(0,Fo)) (138) 
-^F„ + iVF„ = ai^F„_i+R„_i form>l (139) 
where N is the matrix 

r'(ayg(0,Fo)-A) (140) 
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and the function Rm_i(^) depends only on the with k < m: 



(m - 1)/ + ^ 



1 d'^ 



3=0 



z=0 



(141) 



7.5 Formal arguments for thin-film equation 

Consider the particular initial value problem in one space dimension: 

ht+ {h^h^^^)^^0 , h{x,0) ^ 



1 + 



(142) 



This is a special case of ht + {h^hxxx)x — 0. Global existence proofs are available only for 
n > 3.5; numerical solutions suggest finite-time singularity for n — 1 |17| . 

For the problem l|142() and variations of it, the complex region for which existence is expected, 
at least for small t, includes the real a;-axis. For the specific initial value problem, we change 
variables: 

h{x,t) ^ H{l + x^,t), ^ = l + x'^ 

and obtain a nonlinear PDE for H{^,t). A formal asymptotic expansion in powers of t results 
in 



^ 3=0 



(143) 



where P2j are homogenous polynomials of order 2j. With appropriate changes of variables, we 
expect the regularity theorem 5 to be adaptable to prove short term existence for a complex ^ 
sector that includes (l,oo) (i.e., x € R), and show further the validity of (|143|l for ^ 3> t^^^ in 
this sector. 

Asymptotics ()143|l fails when ^ = 0{t^^'). Introducing scaled variables, 
V = it-^'\ T = t^l\ H{^{f^, r), t{T)) = C'F{^, r), 
gives a formal solution as an expansion in integer powers of r, 



(144) 



fc=0 



We expect this series to be convergent. The equation of Fq can be integrated once by using 
far-field matching condition to give: 



^0 ^0 



6 



6 



2^0 ^0 ^f^^f 



3t?" 
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n 



0^0 



Hi 
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= 



With the further transformation ^o('7) = 1 + yiv"^^^), the equation for y is in a form to which 
the general theory [14| applies. From the leading order singularity of the ODE, and the expected 
convergence of ()144|) . as for modified Harry-Dym, we expect to show that the thin-film equation 
has singularities at points close to Xs{t) with 1 + x^ ~ rjst^^'^ . 
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Figure 1 : Subdomains of T) 
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Figure 3: Domain £ in the x-plane 
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